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ABSTRACT 

Taking advantage of the ultra-deep near-infrared imaging obtained with the HST 
on the Hubble Ultra Deep Field, we detect and explore for the first time the properties 
of the stellar haloes of two Milky- Way like galaxies at z~l. We find that the structural 
properties of those haloes (size and shape) are similar to the ones found in the local 
universe. However, these high-z stellar haloes are ~3 magnitudes brighter and present 
bluer colors ((g-r)<0.4 mag) than their local counterparts. The stellar populations 
of z^l stellar haloes are compatible with having ages <1 Gyr. This implies that the 
stars in those haloes were formed basically at l<z<2. This result matches very well 
the theoretical predictions that locate most of the formation of the stellar haloes at 
those early epochs. We note, however, that high-z stellar haloes are slightly (0.5-1 
mag) brighter than what a simple passive evolution of their stars would require to 
evolve them into the local haloes. That suggests that some accretion/star formation 
is still active at z~l. 

Key words: Galaxies: Evolution, Galaxies: Formation, Galaxies: Spiral, Galaxies: 
Structure, Galaxies: Photometry, Galaxies: Stellar Halos 



1 INTRODUCTION 

A natural prediction of the hierarchical galaxy formation 
scenario is the ubiquitous presence of stellar haloes sur- 
rounding the galaxies (e.g. Eggen et al. 1962; Searle & Zinn 
1978; Steinmetz & Muller 1995; Bekki & Chiba 2001; Sam- 
land & Gerhard 2003). These extended and diffuse stellar 
components will be formed from the debris of disrupted 
satellites accreted along the cosmic time (e.g. Brook et al. 
2003; Bullock & Johnston 2005; Abadi et al. 2006). Recent 
cosmological simulations predict that the amount of stellar 
mass in these halos should be ~1O 8 M -1O 9 M Q for Milky- 
Way (MW) like objects. If these simulations are correct, 
most of the stellar halo mass would be assembled before 
z~l and we would expect since then a simple passive evo- 
lution of the stellar populations of these haloes towards the 
present (e.g. Cooper et al. 2010; Font et al. 2011). 

From the observational point of view, the observation 
of the stellar haloes is a major challenge. For example, Mor- 
rison (1993) estimates the surface brightness of the Galac- 
tic halo at the solar radius to be V~27.7 mag/arcsec 2 . In 
addition to the well-characterized Milky Way halo (see a 
recent review by Helmi 2008), deep surveys of M31 (e.g. 
Ferguson et al. 2002; Irwin et al. 2005; Kalirai et al. 2006; 
Ibata et al. 2007; McConnachie et al. 2009) have revealed an 
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extensive halo (to ~150 kpc) with abundant substructure. 
M33 seems also to have a halo similar to the Milky Way 
and M31, despite its smaller total mass (McConnachie et 
al. 2006). Beyond the Local Group, there is growing obser- 
vational evidence showing that stellar haloes are ubiquitous 
and diverse. These studies have been conducted in nearby 
galaxies by using resolved star counts (e.g. Mouhcine et al. 
2005, 2007, 2010; Ibata et al. 2009; Radburn-Smith et al. 
2011) or by extremely deep observations able to detect low 
surface brightness features (e.g. Sackett et al. 1994; Shang 
et al. 1998; Martmez-Delgado et al. 2008, 2009; Jablonka 
et al. 2010; Bakos & Trujillo 2012). Zibetti et al. (2004; 
see also Bergvall et al. 2010), stacking 1047 galaxies, were 
able to observe the properties of stellar haloes up to z~0.05 
(~200 Mpc in distance). Finally, Zibetti & Ferguson (2004), 
using the Hubble Ultra Deep Field (HUDF; Beckwith et al. 
2006) , were able to detect the stellar halo of a disk galaxy at 
z=0.32. This last observation represents the current farthest 
detection of these faint components of the galaxies. 

The observed number and physical properties of the 
streams in the nearby galaxies are in general agreement with 
the predictions from the ACDM model (e.g. Bell et al. 2008; 
Gilbert et al. 2009, McConnachie et al. 2009; Starkenburg 
et al. 2009; Martmez-Delgado et al. 2010). Even the large 
discrepancies in the properties of the stellar haloes of similar 
mass disk galaxies as the Milky Way and M31 have been ex- 
plained, within the cosmological context, as inherent to the 
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system-to-system scatter in the halo formation histories of 
these objects (e.g. Font et al. 2011). To gain, consequently, 
a deeper understanding of the formation mechanisms of the 
disk galaxies is necessary: (a) to explore systematically the 
stellar halos of nearby galaxies to even fainter surface bright- 
ness magnitudes (V>30 mag/arcsec 2 ) where a plethora of 
substructures is expected to appear and (b) to explore the 
stellar halos of high redshift disk galaxies, where the prop- 
erties of the stellar haloes are caught in earlier episodes of 
assembly. The present paper deals with the second possi- 
bility. In particular, if the stellar halo formation scenarios 
are correct, a cosmic epoch that it is worth exploring is z~l 
where, as mentioned before, the haloes should be already in 
place. 

If observing stellar haloes in the local universe is dif- 
ficult, the observation of stellar haloes at z~l requires the 
deepest ever observations from the space. At z=l, the cos- 
mological dimming is ~3 mag/arcsec 2 and the optical red 
side of stellar emission (which contains most of the stellar 
flux) is only visible using near-infrared data. Fortunately, 
the observations to start exploring the stellar haloes up to 
z~l have been recently taken by using the Wide Field Cam- 
era 3 (WFC3) on board the Hubble Space Telescope (HST) 
imaging the Hubble Ultra Deep Field. In this pilot paper, we 
present for the first time the stellar haloes of two MW-like 
galaxies at z~l. We will show that these stellar haloes, in 
agreement with theoretical expectations, were much brighter 
and bluer than present-day stellar halos of similar mass disk 
galaxies. 

The structure of this paper is as follows. In Section 2 
we describe the data used in this paper and the selection of 
our sample, and in Section 3 we conduct a detailed analy- 
sis of our results. Finally, in Section 4 we perform a short 
discussion and a summary of our results. 

All magnitudes in this paper are given in the AB sys- 
tem unless otherwise stated. Throughout we assume a flat 
-dominated cosmology (f2 m =0.3, £}a=0.7, and Ho=70 km 
s" 1 Mpc" 1 ). 



2 DATA AND SAMPLE SELECTION 

Due to the faintness of the features we aim to explore in this 
paper, the study presented here requires the nowadays state- 
of-the-art imaging in terms of depth and resolution. This 
implies using the deepest optical and near-infrared imag- 
ing ever taken with the HST. The Hubble Ultra Deep Field 
(R.A. =03:32:39.0, Dec = -27:47:29.1 (J2000)) is the point 
on the sky with the best HST data for this project. 



2.1 Data 

In the optical regime, the HUDF has been observed with 
the Wide Field Camera (WFC) of the Advanced Cam- 
era for Surveys (ACS) in four niters: F435W (B), F606W 
(V), F775W (i) and F850LP (z). The total exposure 
time in each band is as follows: 134880s (B), 135320s 
(V), 347110s (i) and 346620s (z), and the AB mag. ze- 
ropoints reached are: 25.673 (B), 26.486 (V), 25.654 (i) 
and 24.862 (z). The individual exposures in each filter 
were combined and drizzled to produce a single image per 



band with the following characteristics: 10500x10500 pix- 
els with a pixel scale of 0.03 arcsec/pixel. These data were 
made public by the HUDF team at the following webpage: 
| htt p : / / archive .stsci.edu/pub /hlsp / udf /acs- wf c / 1 

In the near-infrared the HUDF was observed both 
by the Near Infrared Camera and Multi-Object Spec- 
trometer (NICMOS) camera and more recently by the 
WFC3. The WFC3 observations are deeper and with 
higher resolution than NICMOS, however, the WFC3 
observations do not cover entirely the HUDF field. 
For these reasons, in this paper, we have used both 
near-infrared data-sets. The NICMOS HUDF (Thomp- 
son et al. 2005) data consists on F110W and F160W 
drizzled images of 3500x3500 pixels with a pixel 
scale of 0.09 arcsec/pixel. The zeropoints are: 23.410 
( F110W) and 23.22 (F160 W). The images are available 
at |http: / /archive.stsc i. edu/pub/hlsp/udf/ nicmos- treasury /] 
The WFC3 HUDF imaging used here (Bouwens et al. 2011) 
only covers the very central (4.7 arcmin 2 ) part of the HUDF 
but to an extraordinary depth. The zeropoints are: 26.27 
(F105W), 26.25 (F125W) and 25.96 (F160W). The com- 
bined individual exposures in these filters were drizzled and 
have approximately 3000x3000 pixels with a pixel scale 
of 0.06 arcsec/pixel. The data is publicly aavailable at: 
http:/ /archive. stsci.edu/prepds/hudf09/ 

2.2 Sample Selection 

The selection of our galaxies was conducted using the Rain- 
bow Cosmological DatabasfQ published by Perez-Gonzalez 
et al. 2008 (see also Barro et al. 2011a,b). This database is a 
vast compilation of photometric and spectroscopic data for 
several of the deepest cosmological fields, such as GOODS- 
North and South, COSMOS, or the Extended Groth Strip, 
among others. Using all the available photometry, they have 
build spectral energy distributions (SEDs) covering the elec- 
tromagnetic spectrum from the X-ray to the radio wave- 
lengths. Analyzing these SEDs, they have derived very ro- 
bust photometric redshifts and accurate estimates of stellar 
parameters (such as the stellar mass, the UV- and IR-based 
Star Formation Rates (SFRs), the stellar population age, etc 

We have selected within the Rainbow database those 
objects with spectroscopic redshift within 0.8<z<1.2 and 
M*>5xlO 9 M (Kroupa 2001 IMF). This provides us with 
15 galaxies. We have used the ACS z-band imaging to vi- 
sually explore these objects. Two of these galaxies are el- 
lipticals and are not considered on what follows. From the 
remaining 13 galaxies, 9 objects are irregular or appear 
distorted by interactions. There are 4 bona fide spirals: 
UDF3372, UDF4438, UDF5417 and UDF9444. From those, 
we have finally selected only the two disk galaxies with the 
lower inclinations, to mitigate potential effects of the dust on 
the surface brightness profiles we have obtained. The char- 
acteristics of these galaxies are shown in Table 1. We have 
also shown these objects (in z and H bands) in Fig. 1. Unfor- 
tunately, UDF3372 is beyond the area covered by the HUDF 
WFC3 pointing, so we have had to rely our analysis on the 
near-infrared using the NICMOS data. 

1 https://rainbowx.fis.ucm.es/ 
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Figure 1. Our sample of Milky- Way-like galaxies at z~l: UDF5417 and UDF3372. Upper panels: ACS z-band (~B-band restframe) of 
the galaxies in our sample. Bottom panels: WFC3 F160W band (UDF5417) and NICMOS F160W band (UDF3372) (~I-band restframe) 
of our galaxies. Listed on each figure is the galaxy name, its stellar mass, and its spectroscopic redshift. The solid line indicates 1 arcsec 
angular size. 



Table 1. General properties of the selected galaxies 



Name 


R.A. 


Deck 


Redshift 


M abs Stellar Mass 




(J2000) 


(J2000) 




(B-band) (xlO 10 M Q ) 


UDF3372 


03 32 42.3 


-27 47 46 


0.996 


-21.7 2.4 


UDF5417 


03 32 39.9 


-27 47 15 


1.095 


-21.6 3.3 



2.3 Profile Extraction 

We have extracted the surface brightness profiles in all the 
bands available for each of the galaxies. In the case of 
UDF5417 we have obtained both WFC3 and NICMOS pro- 
files, but we only show the deepest profiles corresponding 
to WFC3. The technique used to obtain the surface bright- 
ness profiles is fully explained in Bakos & Trujillo (2012). 
Summarizing, we extract radial surface brightness profiles 
on masked images in order to avoid contamination on our 
light profiles. We apply conservative masking onto sources 
which clearly do not belong to the galaxy, like foreground 
stars, background galaxies, etc. These sources are extracted 
by SExtractor (Bertin & Arnouts 1996). We use some SEx- 
tractor parameters such as the measured flux, elongation, 
and similar to determine the shape and size of these mask 
regions. 

In order to extract radial surface brightness profiles rep- 



resentative of the most external part of the galaxy, we need 
to get characteristic values of the ellipticity and position 
angle of this region. We have done this by computing the 
second-order moment of the light distribution of the galaxy 
using the z-band image. The second-order moment is di- 
rectly related to the position angle, the semi-major (A) and 
semi-minor (B) axis lengths. We fix this ellipticity and po- 
sition angle for all elliptical apertures. In order to increase 
the signal-to-noise counteracting its natural decrement due 
to the lower surface brightness in the outer regime, we used a 
logarithmic radial sampling with steps of 0.03. In each aper- 
ture we estimate the galaxy flux by the 3-cr rejected mean 
of the pixel values of that aperture. This helps to minimize 
the effect of morphological features like a spiral arm crossing 
the aperture. 

Although our images are sky substracted, our light pro- 
files, however, can still be contaminated by some residual 
local sky background. This local background is estimated 
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by using equally spaced apertures. We obtain the number 
count profiles of the galaxy up to very far distances, and we 
chose a large aperture where the profiles become flat, beyond 
the identifiable profiles of the galaxy. The 3-a rejected mean 
of the fluxes inside this aperture gives a robust estimate of 
the residual local sky background (for further information, 
see Pohlen & Trujillo 2006). The error of the background 
determination is key to decide the surface brightness level 
down to which we trust our profiles. Following a conservative 
approach,this is placed where the profiles obtained by either 
over- or undersubtracting the sky measurement by ±<r start 
to deviate with more than 0.2 mag from the original profiles. 

The observed profiles are shown in Fig. [21 The profiles 
have been placed up and down by some magnitudes to fa- 
cilitate their comparison among the different bands. 



3 ANALYSIS 

3.1 The shape of the surface brightness profiles 

The two most remarkable features shown in the surface 
brightness profiles of Fig. [2] are the presence of a break at 
around 1 arcsec (i.e. ~8 kpc at z—1), and some extra light 
above the outer disk expectations starting at R~2 arcsec. 
These two features are observed in all the bands and in the 
two galaxies. This strongly indicates that these features are 
real and are not connected to sky subtraction problems. We 
also note that they are observed in the NICMOS data (which 
are shallower than the WFC3) pointing out that (although 
the profiles are noisier) these data are still deep enough to 
detect the excess of light in the outer regions of the galaxies. 
We have checked that this in fact the case using the NIC- 
MOS F160W filter of UDF5417 and comparing this to the 
WFC3 F160W of the same object. The excess of light is also 
apparent in the NICMOS profile. The position of the break 
is pretty stable (at ~1 arcsec), independently of the differ- 
ent spatial resolution of the different bands used. We note, 
however, that in the very inner region of the UDF3372, the 
presence of a bulge could be masked by the worst resolution 
of the NICMOS imaging compared to the WFC3 data. In 
short, the different resolutions along the various bands only 
show a visible effect in the most inner region of our galaxies. 

We work in this paper under the hypothesis that the 
excess of light at R>2 arcsec is due to the presence of stellar 
haloes in these objects. In order to further investigate this, 
in Fig. [3] we show the result of stacking the images from 
the three WFC3 bands of UDF5417. This stacked image has 
been filtered with a median filter of 5 pixels width to enhance 
the lowest surface brightness features of the galaxy. The re- 
gion corresponding to the excess of light above the outer disk 
expectations observed in the profiles (i.e. 16<R<30 kpc) are 
enclosed by two orange dashed circles. We can ascertain that 
this region is well beyond the main body (disk) of the galaxy. 
We also observe that the distribution of light in that region 
is not perfectly symmetric. This absence of symmetry in a 
region beyond the disk could be an indication of ongoing 
accretion. 

In the local Universe, deep surface brightness pro- 
files (10 magnitudes in range) of spiral galaxies, as the 
ones shown here for z=l objects, also present the excess 
of light beyond the outer disk (see e.g. Bakos & Trujillo 



2012). In present-day spirals, the surface brightness (fi r >28 
mag/arcsec 2 ), colors (g-r>0.6) as well as the radial distances 
(R>15 kpc) of these light excesses correspond with the ex- 
pected characteristics of stellar haloes. On what follows, we 
explore the restframe properties (surface brightness and col- 
ors) of our disk galaxies at z~l. 

3.2 Restframe surface brightness and color 
profiles 

Using the observed surface brightness profiles, we have built 
the SEDs at different radial distances. We have interpolated 
among the different observed filters to obtain the flux in 
the equivalent SDSS restframe bands: u, g and r. After do- 
ing this, we have corrected the restframe surface brightness 
profiles by adding the equivalent magnitude to the cosmo- 
logical dimming at those redshifts. We show the restframe 
surface brightness profiles of our galaxies in Fig. [4] The rest- 
frame surface brightness profiles are used to build the color 
profiles too. For that reason, we show the restframe profiles 
only down to the radial distance where the errors on the 
surface brightness distributions are not larger than 0.2 mag. 

From that figure, it follows immediately that the sur- 
face brightness of our high-z galaxies are brighter (~2 
mag/arcsec 2 ) than those shown in locally equivalent stel- 
lar mass galaxies. This is expected due to the large star 
formation activity and younger stellar population of those 
objects ~8 Gyr ago. Fortunately for us, not only the disk of 
the galaxies are brighter but also the stellar halos. This is 
the reason why we can observe them even considering that 
the surface brightness of these galaxies are severely affected 
by the cosmological surface brightness dimming at z~l (~3 
mag/arcsec 2 ). The radial color profiles of our galaxies are 
also shown in Fig. [4] These color profiles show the U-shape 
(associated to the break in the stellar disk) which has been 
observed in the nearby universe (Bakos et al. 2008) as well 
as in z<l disk galaxies (Azzollini et al. 2008a). In agree- 
ment with the fact that high-z galaxies have younger stellar 
populations at z=l than now, the observed colors of our 
disk galaxies are bluer than in local objects. For instance, 
whereas in the local universe, the surface brightness profiles 
usually have a (g-r)=0.47±0.02 mag at the break radius, our 
galaxies present (g-r)<0.3 mag at this feature. 

In Table 2, we provide some structural values of the disk 
component of the galaxies. In particular, we present the ra- 
dial position of the break and the r-band restframe surface 
brightness (after correcting by the cosmological dimming) 
at this position. Both the values that we get for the sur- 
face brightness of the breaks (^i r ,6ieafe~19.9 mag/arcsec 2 ) 
and their radial positions (Rtreafc~7.7 kpc) are typical of 
the population of disk galaxies with similar absolute magni- 
tudes and stellar mass at those redshifts (e.g. Azzollini et al 
2008b). In the next section we concentrate on the properties 
of the stellar haloes. 



3.3 Stellar halo properties 

Following a similar analysis to the one presented in Bakos & 
Trujillo (2012) for local galaxies, we characterize the struc- 
tural properties of the stellar haloes at high-z by decompos- 
ing the r-band surface brightness restframe profile in three 
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Figure 2. Observed profiles of the UDF5417 (ACS and WFC3) and UDF3372 (ACS and NICMOS) galaxies. The breaks on the disks 
of the galaxies at ~1 arcsec (i.e. ~8 kpc) are observed in all the bands. Beyond R=2 arcsec, in all the cases, there is an excess of light 
over the outer disk (1<R<2 arcsec) light exponential decrement. We identify this excess as the growing effect of the stellar halo in the 
surface brightness profiles of the galaxies. 



Table 2. Structural properties of the selected galaxies 



Name 


^break 








< Me.haio > 


Lhalo/Ltotal 


(g-r)haio 




arcsec 


kpc 


mag/arcsec 2 


kpc 


mag/arcsec 2 


r-band 


mag 


UDF3372 


0.95 


7.6 


19.9 


17 


24.9±0.1 


0.068±0.010 


0.35±0.1 


UDF5417 


0.94 


7.7 


19.8 


19 


24.4±0.1 


0.041±0.010 


0.40±0.1 



different components: a Sersic r 1//n (1968) bulge, a double- 
exponential disk, and a stellar halo described by a Sersic 
model with n=0.5 (see Fig. [5}. This later election is some- 
how arbitrary but it has been chosen because it represents 
well the shape of the profile in the outer region and also be- 
cause it allows a direct comparison with the analysis of the 
stellar haloes conducted at z~0. 

In Table 2, we provide the effective radii of the haloes 
and its mean surface brightness according to the restframe 
r-band. In the same table we also show the fraction of light 
that is contained in the stellar haloes compared to the total 
light of the galaxies (Lhaio/^totai)- The effective radii of our 
haloes (R e ,ha!o~18 kpc) are within the range of the ones 
found in the local universe. The most important difference 



of the high-z stellar haloes compared to the local ones is 
in connection with their mean effective surface brightness. 
At high-z, the values found < /i e ,haio >~24.6 mag/arcsec 2 
(r-band) are significantly much brighter (> 3 mag/arcsec 2 ) 
than those measured in nearby galaxies. We finally have ob- 
tained the (g-r) color at the stellar halo position by averaging 
the different radial colors we get at 16<R<22 kpc (see Fig. 
[4]). Ideally, we would like to explore these colors farther away 
but the errors are so large that this estimation becomes pro- 
hibitive. For this reason, we concentrate on the above region. 
We should, consequently, take these color estimations of the 
stellar halos with caution as they could be slightly affected 
by the colors of the outer disk population. 
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X (kpc) 

Figure 3. Stacked image of all the WFC3 bands (F105W, F125W and F160W) of the UDF5417 (z=1.09) galaxy. The image has been 
filtered with a median filter of 5 pixels width. These allows us to see the faintest features of the object. The dashed blue circle indicates 
the position of the break on the surface brightness profile of the disk of the galaxy. The region enclosed by the two orange dashed circles 
corresponds to the region of the stellar halo explored in this paper. The visible structure of the halo is not perfectly symmetric suggesting 
some on-going accretion. 



3.4 Comparison with the local sample 

To put our high-z galaxies in context, we have compared the 
properties of their stellar haloes with the ones found in the 
local universe. The result of this comparison is presented in 
Fig. [()] In the left panel, we have the light fraction of the 
stellar haloes (in the r-band) versus the total galaxy stellar 
mass. We have taken the local values from different authors. 
These works are listed in the caption of Fig. [5] We find that 
the amount of light contained in the stellar haloes compared 
to the total light of the galaxies is very similar to the values 
found for local reference disks (~5%). This implies that the 
fading with cosmic time of the stellar haloes should be very 
similar to the one that the whole galaxy experiences. 

In the right panel of Fig. [()] we compare the mean ef- 
fective surface brightness and the (g-r) color of the stellar 
haloes both at high and low-z. As we have been claiming 
through the article, the stellar haloes at high-z are signifi- 
cantly much brighter than today. Also, their colors are bluer 
than in the local universe. Can the high-z stellar haloes in 
the high-z universe evolve into the ones that we observe lo- 
cally? The most simple exercise that we can conduct is to fol- 
low the tracks of SSP with cosmic time and to probe whether 
the passive evolution of the stellar population in the haloes 
can reproduce the local properties. These tracks have been 
constructed using the Vazdekis et al. (2010) models for SSPs 
with two different IMFs (Salpeter 1955 and Kroupa 2001). 
We find that the passive evolution moves the color and the 
surface brightness of the high-z stellar haloes towards the 
observed values in the local universe. However, the passive 
evolution is not enough. In fact, the high-z stellar haloes 
are slightly brighter (~0.5-l mag) than what a simple pas- 



sive evolution would suggest. We explore the consequences 
of this finding in the next section. 



4 DISCUSSION AND CONCLUSIONS 

In the current paradigm, stellar haloes of galaxies have 
a dual nature. The outer halo formed primarily by ac- 
cretion/mergers, whereas the inner halo is formed mainly 
through dissipational collapse (as a result of gas-rich merg- 
ers at early times). Recent cosmological simulations (e.g. 
Cooper et al. 2010; Font et al. 2011) indicate that both phe- 
nomena take place at z>l. If this is the case, one would 
expect that the stellar haloes will evolve since that epoch 
by mostly fading without any significant change in their 
structures. Our present observations allow us to check this 
important prediction of the disk galaxy formation scenario. 

According to our data, the structural properties of our 
high-z stellar haloes seem to be at place at z~l. Both the 
shape and the sizes of these components are very similar 
to the ones found in the local universe. The major differ- 
ence with respect to the local galaxies is that our haloes 
are ~3 magnitudes brighter than their local counterparts. 
They also present bluer colors than the local stellar haloes. 
Both results agree with the idea that high-z stellar haloes 
are significantly much younger than present-day ones. Ac- 
cording to the SSP evolutionary tracks that we have ex- 
plored, the colors of our high-z stellar haloes correspond to 
stellar populations with an age <1 Gyr. That implies that 
the stars in those haloes were formed basically at l<z<2. 
This matches very well the theoretical predictions that move 
most of the formation of the stellar haloes at those epochs. 
We note, however, that our stellar haloes are slightly (0.5-1 
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Figure 4. Upper panels: u, g and r-band equivalent restframe surface brightness profiles obtained by the interpolation of the observed 
band profiles. The restframe surface brightness profiles have been corrected by the cosmological dimming. The profiles are only shown 
until the error bars reach 0.2 mag. Bottom panels: The g-r and u-g restframe color profiles are shown. The characteristic U-shape color 
profile (very clear in the g-r color) is seen. The bluest color is found at the break position of the disk galaxy. The color of the stellar 
haloes are among the reddest that can be found in those galaxies. 



mag) brighter than what a pure passive evolution of their 
stars would require to evolve them into the local haloes. That 
suggests that some accretion/star formation is still active at 
z~l. 

Our sample is so far very small (only 2 objects), but 
this pilot work can be easily extended towards other spiral 
galaxies in the HUDF at lower redshifts and with different 
inclinations. Our two disk galaxies have not been preselected 
to have bright (i.e. detectable) stellar haloes but we admit 
we can not infer, with only two objects, whether this is the 
general behavior of the whole disk galaxy population at z~l. 
It will be interesting to explore in the future whether other 
disk galaxies have prominent stellar haloes as the ones dis- 
cussed in this paper or whether their stellar haloes are fainter 
and their evolution better fit with a pure passive evolution 
of their stellar population down to the present-day universe. 
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